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Abstract. In this paper, we discuss the equilibrium phases and collapse transitions of a lyotropic nematic 
gel immersed in an isotropic solvent. A nematic gel consists of a cross- linked polymer network with rod- like 
molecules embedded in it. Upon decreasing the quality of the solvent, we find that a lyotropic nematic gel 
undergoes a discontinuous volume change accompanied by an isotropic-nematic transition. We also present 
phase diagrams that these systems may exhibit. In particular, we show that coexistence of two isotropic 
phases, of two nematic phases, or of an isotropic and a nematic phase can occur. 
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1 Introduction 

A polymer gel consists of flexible polymer chains that 
are cross-linked to form a three-dimensional network El 
whose volume in response to solvent changes can swell 
by factors as great as 1000 or more and undergo dis- 
continuous changes along a coexistence line terminating 
in a liquid-gas like critical point [Q. Dispersions of long 
rigid rods undergo a transition from an isotropic to a ly- 
otropic nematic phase with increasing rod volume frac- 
| tion HQ. In this paper, we will explore the possibility 
of creating lyotropic nematic phases in polymer gels with 
rigid rods either dispersed in their open spaces or form- 
ing parts of their constituent polymer chains. Our basic 
premise is that decreases in solvent quality contract the 
gel, thereby increasing rod volume fraction and inducing 
a transition to a nematic phase. A lyotropic nematic gel 
has the same spontaneous broken rotational symmetry 
and macroscopic elasticproperties as a thermotropic ne- 
matic or elastomer J§H(t]|| . As a result, it will exhibit the 
characteristic "soft" elasticity |8|P, ^0| , pd| of a nematic gel 
whereby the shear modulus for shears in planes containing 
the anisotropy axis vanishes and stresses up to critical dis- 
tortions vanish for certain elongational and compressional 
strains. This "soft" elasticity along with unusually large 
coefficients of thermal expansion suggest interesting tech- 
nological applications for thermotropic nematic or other 
liquid crystalline gels including artificial muscles |I^ , |l3| |, 
actuators |Q, and electro-optical devices |jl5). Lyotropic 
nematic gels might find use in similar applications. It is 
also imaginable that nature has found ways to use biologi- 
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cal versions of these materials to control elastic properties 
of cells. 

A theory for lyotropic nematic gels should be sophis- 
ticated enough to treat both the volume-collapse transi- 
tion [ ^6|Jl7|l of an isotropic gel and the development of 
nematic order. Tanaka |2| explained the former phenom- 
ena using Flory's theory for gels JT], p~8|JT9| , which combines 
the Flory-Huggins theory of polymer solutions with the 
Flory theory for the entropic elasticity of isotropic elas- 
tomers. Brochard po|| considered gels in thermotropic liq- 
uid crystal solvents and concluded that a gel crosslinked in 
the isotropic phase of the solvent collapses when the sol- 
vent undergoes a transition from the isotropic to the ne- 
matic phase. Warner and Wang modified Flory's the- 
ory to treat anisotropic gels using the neoclassical theory 
[pd[ for the elasticity of anisotropic rather than isotropic 
elastomers. They then used the phenomenological Maier- 
Saupe theory p2j to describe the development of nematic 
order. Coupling between elasticity and nematic order is 
provided by the dependence of anisotropic step-lengths 
in the neoclassical elasticity theory on the nematic order 
parameter. They found a rich phase diagram, with tem- 
perature as a control parameter, in which coexistence of 
nematic and isotropic phases and a discontinuous uptake 
of solvent are possible. In this paper, we follow Warner 
and Wang's treatment of gel elasticity, but we use the On- 
sager theory ||, which provides a quantitatively accurate 
description for the isotropic-to-nematic (I-N) transition in 
dispersions of hard rods [I| , rather than Maier-Saupe the- 
ory to describe the development of nematic order. Our 
theory provides a more accurate treatment of the strongly 
first-order I-N transition of dispersions of hard-rods than 
does one based on the phenomenological Maier-Saupe the- 
ory, which is best suited for a description of the vicinity 
of the transition in thermotropic systems. 
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(a) 



(b) 



Fig. 1. Representation of the nematic gels referred to as (a) 
model 1 and (b) model 2. In both cases, rods of diameter a 
and length b are embedded in a crosslinked polymer network. 
In model 1, the rods are part of the network which consists of 
crosslinked polymer chains composed of rigid rods and flexible 
spacers. The thickness of the spacers is equal to the diameter 
of the rods, and the rods are much smaller than the average 
distance between crosslinks. In model 2, the independent rigid 
rods are dispersed in the free spaces of a crosslinked flexible 
polymer network. There is no particular constraint on the dis- 
tance between crosslinks and on the size of the chains which 
constitute the network. 



We consider two models depicted schematically in Fig. 
Ill; model 1 in which the polymer chains constituting the 
gel consist of rigid rods connected by flexible segments 
and model 2 in which rigid rods are dispersed in the fluid 
volume between polymer segments. Though we introduce 
certain not completely controlled approximations to be 
discussed below, we believe that these models provide a 
realistic semi-quantitative description of lyotropic nematic 
gels. In both models, we find as expected a first order I- 
N transition accompanied by a discontinuous change in 
volume. As in the systems considered by Brochard p(| , 
the development of nematic order is accompanied by a 
decrease in the gel volume. In model 2, we also find lines 
along which either two isotropic or two nematic phases of 
differing densities coexist. These lines terminate in mean- 
field liquid-gas-likc critical points, which are expected to 
retain their mean-field character even when fluctuations 
beyond mean-field theory are included because of effec- 
tive long-range forces induced by a non-vanishing shear 
modulus §. 

This paper is organized as follows: In the next section 
(Sec. ||), we develop the model free energy for both models 
1 and 2 including the orientational free energy of rods 
described by the Onsager theory. In Sec. [|, we construct 
the phase diagrams for models 1 and 2. We conclude with 
a brief summary. 



2 Free energy of nematic gels 

In this section, we will derive the free energy for nematic 
gels. This free energy includes entropic and enthalpic con- 



tributions that are present in uncrosslinked polymer solu- 
tions, clastic contributions arising from crosslinking, and 
contributions arising from orienting rigid rods, and it de- 
pends the volume fractions of polymer, solvent, and rods. 
We denote by cf> the total volume fraction of the network, 
which contains flexible polymers of volume fraction (f>p 
and rods of volume fraction 4>r, so that 



OR- 



(1) 



The solvent volume fraction is thus 1 — <fi. Under the as- 
sumption that the rods are confined within the network 
as it deforms, the ratio of the concentration of the rods to 
that of the network must be a constant, which we denote 

by 



(2) 



In model 2, we have not considered the possibility of an 
equilibrium of the mobile rods between the inside and the 
outside of the gel as done in Ref. |2(| for flexible rods 
chemically identical to the polymer chains of the network. 
Instead by using Eq. with a fixed /3, we have assumed 
that a constant number of rods stay trapped within the gel 
as the gel changes its volume. It is possible that the rods 
stay confined in the polymer network either for kinetic 
reasons or because of entanglements within the network. 

To describe nematic gels immersed in a solvent, we 
assume that the total free energy can be written as a sum 
of three contributions 



J~ — J~ ri 



•F °l ~i~ *F rod 



(3) 



where !F m ix is the free energy of mixing of the gel with its 
solvent, J- e i is the elastic free energy of the gel network, 
and T ro d is the free energy of the rods. Within the Flory- 
Huggins theory, T m i X may be written as ||l|,p|,fjl| 



3~> n 



Vk B T 



fn 



where V is the volume of the gel and v c 
a lattice site, which is assumed to be v c 



(4) 

is the volume of 
= a 3 , where a is 



the monomer diameter, 
site given by 



f m ix is the free energy per lattice 



fmix = (l-0)ln(l-0)+ X <Ml-0)+Xfl.<Ml-0) 
+ Xp<Pr<Pp, (5) 

where \i XR> an d Xp are the Flory-Huggins parameters 
that characterize, respectively, the polymer-solvent, the 
rod-solvent, and the rod-polymer interactions at the level 
of the second virial approximation. The first term in Eq. (||) 
represents the entropy of the solvent. Note that the trans- 
lational entropies of the polymers, 0pln0p, and rods, 
ipRlnipR, are absent in Eq. dq) because the chains, which 
include the rods, are attached to the network fig] . In or- 
der for the Onsager theory to be valid, the solvent quality 
must be good for the rods so that their mutual interactions 
are well described by the excluded volume interaction. For 
simplicity, we assume that xr = Xp = corresponding to 
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an athermal solvent for the rods and absence of a direct 
rod-polymer interactions. Thus, the single parameter Xi 
which measures the strength of interaction between poly- 
mer and solvent, controls the solvent quality. It is expected 
to vary with temperature as 0/2T, where is the theta- 
temperature separating good-solvent (T > 6*) from poor- 
solvent (T < 0) behavior. 

The elastic free energy T e i introduced by Warner et al. 
]Tl[ describes the energy change associated with a defor- 
mation of a nematic elastomer characterized by homoge- 
nous strains parallel A|| and perpendicular Aj_ to the direc- 
tor (along which the rods are aligned). It may be written 
as 



the swelling ratio by 



T A = h B TN c 



A2 e 



(0) 



where N c is the total number of network strands, £q is the 
step length of the polymer in the isotropic phase in which 
we assume the system was cross-linked, and £m and £± are, 
respectively, the step lengths parallel to and perpendicu- 
lar to the director in the nematic phase. Strains A arise 
spontaneously due to the anisotropy of the network, and 
they are measured with respect to an isotropic state char- 
acterized by £q. Note that in and £± are in general func- 
tions of the de Gennes-Maier-Saupe |^2| order parameter 
Qij , which describes the orientational ordering of the rods, 
and they reduce to £q if Qij = 0. Equation (||) is based 
on the following assumptions: (i) all chains are sufficiently 
long that their end-to-end vectors may be described by an 
anisotropic Gaussian distribution, and (ii) this distribu- 
tion transforms affinely under macroscopic deformation. 
Equation (|^) has been used to account for many proper- 
ties of nematic elastomers 0,^1 23 . 

Equation (^) differs from the Flory expression for the 
elastic free energy of a gel QJl^l , which contains an addi- 
tional strain-dependent term of the form In (a|A_l). While 

Tanaka et al. Q and Matsuyama et al. (24j have used the 
Flory theory with this term included to explain the col- 
lapse of isotropic gels, we believe (and later check) that 
this term does not play a significant role in nematic gels. 
Therefore, we do not include this term in the discussion 
that follows [ ^TlJl^ , ^5| . It is, nevertheless, interesting to 
note that this term is expected to be absent in the so-called 
phantom model of rubber elasticity |2{|. This model es- 
sentially ignores the effect of entanglements of the chains 
and thermal fluctuations of the junctions points and is 
believed to be valid in the dilute regime. In more refined 
models, however, where these effects are properly taking 
into account, there may be additional strain dependent 
terms contributing to Eq. (^|) p6| . 

Unlike thermotropic nematic elastomers, which are es- 
sentially incompressible, lyotropic nematic gels can un- 
dergo large volume changes because of their ability to 
take up solvent [p7| . We measure volume changes by the 
swelling ratio V/Vo — </>o/</>, where Vq is the network vol- 
ume (and 4>o its volume fraction) in the reference state, 
and V is the network volume (and <f> its volume fraction) 
in any particular state of gel. The strains are related to 



V 



An A i — — — 



V <t> 



(7) 



This relation allows us to express the elastic free energy 
Eq. (|^) in terms of A|| only. Minimizing Eq. (JsJ) with re- 
spect to An, we obtain 



*ll = 



"Mil 



and 
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where v c = N c v c /Vq. Note that Eq. (||) is slightly more 
general than the one used by Warner et al. pl[ , where 
the particular case oi 4>q = 1, i.e. & dry gel, is considered. 
As pointed out in Ref. [Q, in order for an isotropic gel to 
undergo collapse it is crucial that </>o <C 1, and we expect 
a similar conclusion to apply to nematic gels. 

The free energy T ro d describes the physics of the rods 
embedded in the gel. Unlike previous treatments of ne- 
matic elastomers and gels, which focus mostly on ther- 
motropic systems, our current treatment focusses on ly- 
otropic gels. We employ Onsager theory |3| to describe 
the isotropic-nematic transition that occurs upon increas- 
ing rod volume fraction via osmotic compression of the 
gel. The free energy obtained by Onsager is H 



= N R J d 2 u*(u)hi 



Ve 



ab 2 N R 
V~ 



J d 2 udV<F(u)#-(u')|u x u'|, (10) 



where u is the unit vector specifying rod direction, N R 
is the number of rods in the network with length b and 
diameter a and ^(u) is the orientational distribution func- 
tion. The first term in Eq. flio| ) is the entropy of the rods, 
and the second takes into account the excluded-volume 
interaction between the rods at the level of the second 
virial approximation, which becomes exact in the limit of 
infinitely long rods. Equation (jlTJ) is the Onsager free en- 
ergy for rods in an isotropic solvent, and not for rods either 
attached to or dispersed in a network. It clearly overesti- 
mates both rotational and translational entropy of rods 
since the network will introduce steric constraints whose 
effects are difficult to estimate. Using the Onsager free en- 
ergy for the rods rather than one that takes full account 
of constraints imposed by the network is the principal ap- 
proximation of our theory. As discussed below, we will 
treat the rod translational entropy differently in models 1 
and 2. 

Even in the absence of the coupling between the strains 
and nematic order, minimizing the free energy J- ro d yields 
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an intractable integral equation for !?(u). We, therefore, 
follow Onsager Ppq] and introduce a trial distribution 
function of the form 



*(u) 



acosh(au • c) 
4-7T sinh a 



(11) 



where a is a variational parameter and c is an arbitrary 
unit vector. Thus, the two integrals in Eq. (|l(]) can be 
performed analytically: 



a(a) 



d 2 uV(u) In [AttW(u) 



In (acotha) — 1 



arctan (sinh a) 



p(a) = 



sinh a 

- J d 2 udV ^(11)^(11') |u x u'| 

2/ 2 (2a) 



sinh a 



(12) 



(13) 



where l2{x) is a modified Bessel function. 

Using these results and Eq. we can calculate the 
rod free energy per site. In model 1, the rods are locked 
to the network, and there is no independent translational 
entropy associated with them, although there is rod trans- 
lational entropy in f e i arising from fluctuations of the 
polymer-rod segments constituting the network. In model 
2, rods are not attached to the network, and we assume 
for simplicity that they are free to diffuse throughout the 
volume of the gel not occupied by polymer segments and 
that they therefore have a translational entropy per site 
of 4>r In <f)R. In real systems, rods may in be sterically con- 
fined to finite regions of the gel in which case, our model 
overestimates the rod translational entropy. The free en- 
ergies fl od and f 2 od are thus, respectively, 



fr* 



rod 



•> rodVc 

k R TV 



^Pa(a) + -<t?(i 2 p(a) 
irb 7r 



4a 



frod = frod+ -T</>/?ln< 

irb 



(14) 
(15) 



up to an irrelevant constant in the chemical potential. 
Note that the term proportional to <p In <f> represents the 
translational entropy of the rods. The order parameter of 
the I-N transition is defined by the tensor [E9| 



Qij = J d 2 uV(u) 



j 3 ^ 



(16) 



This tensor can be diagonalized to yield 

Q = Diag [-5/3, -S/3, 25/3] , 
where S is the scalar order parameter given by 



To complete our description of nematic gels, we must 
specify the relation between the anisotropy of the polymer 
chains and the nematic order. We define a step length ten- 
sor by iij — (RiRj)/N£ , which conveniently represents 
the anisotropy of the chains, where R is the end-to-end 
vector of the chains between two crosslinks J3(J and N is 
the total number of elements in these chains. The eigen- 
values of are i\\ and £j_, which enter into the elastic free 
energy Eq. (||). In general, £y is some unknown function 
of the nematic order parameter , which may be model 
dependent. The simplest form of this function allowed by 
symmetry is 



1 = ^4" 



(18) 



with an unknown coupling parameter ?/, which depends 
on the network structure and which will be different for 
models 1 and 2. Such a linear relationship is well verified 
experimentally for nematic elastomers pi] . 

In the case of model 1, we assume that there are jN 
rods of length b and volume 7ra 2 6/4, and (1 — 7) TV spacers 
of size a and volume a 3 , which are distributed randomly in 
a polymer chain, forming a main-chain liquid crystalline 
polymer [ S3 . Using the freely-jointed chain (FJC) model 
@||, we find that 



V = Vi 



67(6/0) 



1 - 7 + 7(Va) 2 



(19) 



where 7 is related to (3 for model 1 in the following way 



= 



7T&7 



4a (1 — 7) + 7T&7 



(20) 



Note that Eq. (|l9| ) is exact for the FJC model and is valid 
for any form of the distribution function tf'(u). 

For model 2, on the other hand, the FJC result is 
not expected to apply since the rods are not physically 
connected to the network chains. However, it is reason- 
able that for a sufficiently low concentration of the rods, 
Eq. (|lj) should still hold with an unknown coefficient 77, 
which we take arbitrarily to be 77 = = 6 (the maximum 
value of r] in model 1). Using Eqs. (O) and (|i~8|), we have 



£ \\( a ) 1 , V 0/ X 

—„ — = 1 + iM"). 



l-|5(a) 



Il(cQ 



(21) 
(22) 



In both models S > in equilibrium in the nematic phase, 



S(a) = - /d 2 u!/?(u) [3(u • c) 2 — l] ^11 > £j-> an d l ne adopts a prolate shape. Generaliza- 

2j 1 tion to allow an oblate phase are not difficult but will b( 



1 3 ( 1 1 
1 coth a 

a \ a 



(17) 



Note that S(a) vanishes in the isotropic phase when a = 
and that S(a) > in the nematic phase when a > 0. 



tion to allow an oblate phase are not difficult but will be 
left for another study. Our final free energy per site for 
model 1 is f m ix + fel + flod with the parameter 77 given by 
Eq. ( pj| ) while that for model 2 is f m ix + fel + frod w ith 
rj = 7/2 =6. We analyze the phase diagrams arising from 
these free energies in the next section. 
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3 Equilibrium swelling of a nematic gel 

In this section, we discuss the equilibrium properties of a 
nematic gel. The orientational order of the nematic gel 
is controlled by the parameter a, which is determined 
by the condition that the total free energy per site / = 
/rote + fel + frod is a minimum at a — a*{4>). The con- 
straint <f> < 1 implies that a*(<p) < a max , where a max 
is the value of a* obtained for <j) — 1. At a — a max 
the maximum possible value of scalar order parameter S 
is reached. Note that the orientational order is governed 
only by f rod and f et and not by f mix . The elastic free 
energy couples to the nematic order through £\\(a) and 
£±(a) as given in Eqs. ( pl| ) and (p2|). As pointed out in 
Ref. [ pl[ , although the contribution from f e i is not neg- 
ligible, the presence of the cross-linked network does not 
qualitatively alter the nematic properties of the rods. In 
particular, if the volume fraction <fi is small, / has only one 
minimum at a = 0, which corresponds to the isotropic 
state. If > 4>i, another minimum appears at positive 
a, which corresponds to the nematic state. If <fi > 02, 
the minimum at a = disappears and there is only one 
minimum at a positive a. Therefore, the isotropic state is 
only possible for < (f> < 02 and the nematic state for 
01 < < 1- These features are also found in Onsager 
theory for lyotropic liquid crystals. However, in nematic 
gels, the equilibrium phases are more interesting due to 
the coupling of elastic strains and nematic order. 

The chemical potential fi and the osmotic pressure 77 
are two quantities characterizing a nematic gel. They are 
respectively, given by 



relax, and its concentration and volume are thus fixed. 



H = df/d<f>\ a=a * w , 
dT _ _k B T 
~dV ~ v7 



n 



f-<t> 



df 



(23) 
(24) 



a=a* ($>) 



In the nematic phase, their dependance on includes con- 
tribution from the dependance of a*(0) on 0. The chem- 
ical potential governs the uptake of the solvent molecules 
and the osmotic pressure the mechanical stability of the 
gel. If 77 > 0, the gel is stable with the surrounding sol- 
vent. If 77 < the gel becomes unstable in the absence 
of an external applied stress, and consequently it shrinks 
and separates itself from its fluid until it reaches 77 = 0. 
Thus, the equilibrium swelling of the gel follows the vol- 
ume curve 4>*(x) defined by the condition 77 = 0, which is 
equivalent to minimizing the total free energy T with re- 
spect to 0. Note that in analogy with the phase diagram 
of the liquid-gas model, it is interesting to consider iso- 
bars 0*(x; 77o) describing a gel to which a constant pres- 
sure 77o is applied (l). In the following, we focus mainly 
on the situation with no applied stress 77o — 0, and we 
will not consider situations in which 77o < 0, which can 
only be reached by attaching the outer surfaces of the 
gel to movable walls that can be subjected to tension. In 
order to reach an equilibrium state of the gel, the exter- 
nal condition, such as the temperature, must change very 
slowly. On the other hand, if the temperature, and thus 
X, changes rather rapidly, the gel does not have time to 



In this case, macroscopic phase separation 16 may take 
place if 4>a < 4> < 4>b, which means that an isotropic (or 
nematic) phase of volume fraction 4>a coexists with an 
isotropic (or nematic) phase of volume fraction 4>b within 
the same gel, where 4>a and 4>b are determined by the 
common tangent construction (i) the osmotic pres- 
sure in the two phases are equal: 77(0^) = 77(0^) and 
(ii) the chemical potential are equal: h{4>a) = m(0b)- We 
emphasize that in Onsager theory, the I-N transition is in- 
dependent of temperature and depends only on the volume 
fraction of the rods. Therefore, in our model for nematic 
gels, x is t ne only parameter that depends on the temper- 
ature (x ~ 1/7 1 ). This is in contrast to previous studies 
where the temperature dependence comes from the 
Maier-Saupe interaction parameter and x is consequently 
assumed to be a constant independent of temperature. We 
present in what follows the phase diagrams for models 1 
and 2 in the coordinates vs. % or 77 vs. x f° r fixed values 
of the parameters v c , 0, b/a, and 0o- 



3.1 Nematic gels with rods linked to the network: 
model 1 

In model 1, the rods are part of the network in which a 
polymer chain consists of rigid rods and flexible spacers. 
Note that the parameters 4>o,i; c , f3,b/ 'a and £q are not all 
independent. This is because we have assumed that the 
deformation of the network is affine in deriving the elastic 
free energy Eq. ([}[). This in turn implies that there is only 
one characteristic length scale of the network, the end-to- 
end separation of its constituent chains, and this length 
must be of the order of the average distance between cross- 
links H|. This condition imposes that v c ~ (a/£ ) 3 N~ 3 / 2 
and 4>o ~ (a/io^N^ 1 / 2 , where N is the total number of 
elements (spacers+rods) in the chains between cross-links 
||. Note that only when TV ^> 1, the average distance 
between cross-links is larger than the rod length, which is 
required to justify a posteriori Eq. (||). For example, for 
b/a = 50, (3 = 0.284 and TV = 100, we find that v c ~ 
7 • 10~ 6 and <fio ~ 5 • 10 -3 . These values are used below to 
construct the phase diagram. 

In Fig. ^ the phase diagram of the gel is given in two 
sets of coordinates: in (a) in terms of <fi vs. x an d O 3 ) i n 
terms of the osmotic pressure 77 vs. x- Thick solid lines 
are isobars of constant osmotic pressure for 77 = 0.005, 
and 77 = 0. The most remarkable feature of these iso- 
bars is the discontinuity in volume fraction (collapse) at 
a particular value of x defined as x* ■ The isobar 77 = is 
called the volume curve, and it delimits the region of ab- 
solute stability of the gel in the absence of applied stress. 
Collapse on the volume curve occurs at x* — 0.94. We em- 
phasize that this collapse is fundamentally different from 
the collapse of an isotropic gel discussed by Tanaka . In 
our case, the collapse is driven by the first-order nematic- 
isotropic transition as can be seen in Fig. [|, which shows 
the discontinuous change of the nematic order parameter 
S along the volume curve. Note that S ~ 1 in the nematic 
phase, which means that the rods inside the gel are highly 



D. Lacoste, A.W.C. Lau and T.C. Lubensky: Phase Transitions in Lyotropic Nematic Gels 



(a) 




n =0.006 



0007 



0.7 0.8 0.9 1 

X 



1.1 



(b) 



0.006 i 




Fig. 2. Phase diagram for model 1 in the coordinates: (a) vol- 
ume fraction of the network <j> vs. Flory-Huggins parameter \ 
and (b) osmotic pressure 77 vs. \- Thick solid lines are iso- 
bars of constant osmotic pressure for 77 = 0.006, 77 = and 
77 = 0.0007. Thin solid lines delimit the isotropic-nematic co- 
existence: the lower line 4>a(x) is the boundary of the isotropic 
phase (I) and the upper line 4>b(x) is the boundary of the ne- 
matic phase (N). The isobar 77 = is called the volume curve, 
and delimits the region of absolute stability of the gel in the 
absence of applied stress. The volume curve is discontinuous 
(collapse) at x* — 0.94, which corresponds to the point d of co- 
ordinate (x* 5 0) in (b). The parameters are b/a = 50, j3 = 0.284, 
v c = 7 ■ 10" 6 and O = 5 ■ 10" 3 . 



ordered. This is consistent with experimental observations 
on nematic elastomers |3l| . As shown in Fig. ||, the devel- 
opment of nematic order in the gel results in anisotropy 
in the strains. 

In Fig. ||a, thin solid lines delimit the isotropic-nematic 
coexistence: the lower line <\>a (x) is the boundary of the 
isotropic phase and the upper line 4>b(x) is the boundary 
of the nematic phase. As mentioned in Sec. |^, if the exter- 
nal condition changes rapidly and the system is quenched 
to a particular point within the coexistence region, macro- 
scopic phase separation will occur. The volume fractions 
of the rods in the two phases are /3(f) a and (3(\>b , respec- 
tively. They are both of the order of a/b, which is smaller 



0.8- 



0.6- 



0.4- 



0.2- 



u 0.5 0.6 0.7 0.8 0.9 1 1.1 1.2 

X 

Fig. 3. Scalar nematic order parameter S as function of the 
Flory-Huggins parameter \. The volume curve is the thick line 
and the thin curve is the nematic branch of the isotropic and 
nematic coexisting phases of Fig. ^[ All the parameters are the 
same as in Fig. ^. 




X=0.6 




X=0.94 



Fig. 4. Schematic picture of the extension of the chains inside 
the gel as characterized by the strains Ay and Aj_, correspond- 
ing to the parameters of Fig. ^. The upper picture corresponds 
to x — 0-6 in the isotropic phase, and the two lower pictures 
correspond to the isotropic and nematic phases which coexist 
when the collapse occurs at x* — 0.94. After the collapse, the 
resultant strains are An ~ 0.72 and Ax ~ 0.17. 
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than the maximum volume fraction of the rods f3. Note 
that the coexistence region widens as x increases. This 
may be attributed to the interaction between the rods 
becoming less repulsive as \ increases. A similar feature 
has been pointed out by Flory [Q and Flory and Warner 
fjsj in a different system consisting of a solution of rods 
in the absence of any elastic network. Note also that a 
qualitatively similar behavior has been predicted theoret- 
ically |2^,^6| and observed experimentally [Q in solutions 
of rods where the interaction between rods is attractive 
and is induced by depletion interaction. 



3.2 Nematic gels with mobile rods: model 2 

In model 2, the rods are mobile but confined within the 
gel. Unlike model 1, the parameters in model 2 are less 
restricted. Using a typical value for the shear modulus of 
an elastomer ~ 1CU 5 J/m 3 [^3) and from Ref. we es- 
timate that Vaj^ 1CP 5 and <j)Q ~ 0.001. As can be seen 
in Figs. H to {?] with different values of v c , (3, and b/a, 
the phase diagrams of model 2 have two common features 
with model 1: the discontinuous collapse and the coexis- 
tence of isotropic-nematic phases. In addition, model 2 has 
two additional interesting features: the coexistence of two 
isotropic phases along a line terminating at a critical point 
and the coexistence of two nematic phases along another 
line terminating at another critical point. 

To see the coexistence of two isotropic phases, we ex- 
pand the osmotic pressure Eq. (|2~i| ) in the isotropic phase 
(a = 0) up to fourth order in 0: 



(a) 



nv c 

k B T 



(25) 



where S = | 



- x(l + 4/? 2 A and v c = vjcp 1 ^. The 
last term £0 in Eq. (p5|), with £ = 4a/?/7r6 comes from the 
translational entropy of the rods proportional to /3</>ln^> in 
Eq. (|l5|), which is specific to model 2. The coexistence of 
two isotropic phases ends at a critical point, the location of 
which is found by requiring the first and second derivative 
of the osmotic pressure with respect to <j) to be zero 77' = 
77" = 0. This critical point located at (6 c ,4>c) only exists in 
fact when 77(<5 C , <fr c ) > 0. This last condition imposes that 
£ > £ e , with £ e = ^Zl^^v,?^. Therefore as £ varies, 
the critical point follows a line which ends at a critical 
endpoint when £ = £ e j|. Since £ e > £ = in model 1, 
the volume curve does not have a critical point in this 
case and the coexistence between two isotropic phases is 
not possible, as found indeed in the previous section. 

Equation (|25| ) was derived by Tanaka for an isotropic 
gel with £ = v c /2(j>Q. As discussed in section 0, the term 
proportional to £ arises from the additional strain depen- 
dent term In ^AyA_i_J in the elastic free energy. Hence, the 

existence of the line of coexistence of two isotropic gels 
terminating at a critical point as observed by Tanaka is 
possible only if this term is included In the absence 
of this term, there could only be a continuous volume 
change induced by the change of quality of the solvent. 




0.01 i 



3.5 



4.5 



5.5 



X 



Fig. 5. Phase diagram for model 2 in the coordinates: (a) 
volume fraction of the network <f> vs. Flory-Huggins parameter 
X and (b) osmotic pressure 77 vs. \- The thick solid line in 
(a) is the volume curve and thin solid lines are coexistence 
lines of isotropic-nematic and isotropic-isotropic phases. In the 
region denoted A in (a), the osmotic pressure of the gel is 
negative. Solid lines in (b) represent first order transitions. 
The parameters are b/a = 8, /3 = 0.628, v c = 5 ■ 10 and 
4>o = 10 -2 . There is a discontinuous volume change (collapse) 
at x* — 5.68, between an isotropic phase (I) of volume fraction 
cj>i SO ~ 1.1 • 10 -4 and a nematic phase (N) of volume fraction 
4>nem = 0.9. The coexistence region of the two isotropic gels 
is denoted 7i + 12 and terminates at a critical point denoted c 
and located at \ c ~ 3.67. At \t — 4.31 there is a triple point 
represented by the dotted line in (a) and by the point t in (b) , 
where two isotropic gels of different volume fractions coexist 
with a nematic gel. 



In contrast, we have found in the case of nematic gels a 
discontinuous volume change even without including the 
controversial strain dependent term. Thus, this confirms 
that the physical origin of the collapse for nematic gels is 
the isotropic-nematic transition. 

Fig. |E| shows the coexistence lines between two isotropic 
gels terminating at a critical point at Xa in the case where 
this critical point is located at Xc < X* ■ The dotted line 
in Fig. o represents a triple point, where the two isotropic 
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A 







2 2.5 3 3.5 4 4.5 5 
1 



Fig. 6. For model 2, volume curve (thick line) as a function 
of the Flory-Huggins parameter \, an d coexistence region of 
isotropic- nematic (thin lines). The parameters are the same as 
in Fig. ^| except for v c which has been increased to 7 ■ 10 -4 . 
Now the collapse is at \* — 4.1. 



gels coexist with a nematic gel. For a larger value of v c , 
we note that either \ c > \* or when there is no critical 
point, and the phase diagram looks very much like Fig. 
as shown in Fig. ^. Thus, the isotropic-nematic coexis- 
tence region has swallowed the region of coexistence of 
two isotropic gels. This means that in this region phase 
separation into two isotropic phases is immediately fol- 
lowed by phase separation into a nematic phase and an 
isotropic phase. Note also that in this case, the collapse 
occurs at a smaller value of x* — 4. This is expected since 
the elasticity of the network favors contraction of the gel. 
For completeness, the phase diagram of Fig. [5] is shown in 
the coordinates osmotic pressure 77 vs. \. 

When the ratio of the volume fraction of the rods and 
that of the network defined as j3 is increased while keep- 
ing the other parameters fixed, the coexistence region of 
isotropic and nematic phases moves to lower volume frac- 
tion cf>. This opens up the possibility of coexistence of two 
nematic phases at higher volume fraction with different 
volume fractions and different order parameters. This is 
indeed the case as shown in Fig. I for (3 = 0.78 and 
b/a — 100. The coexistence region of the two nematic 
phases ends at a critical point located at x'c — 3-1 and 
<f>' c ~ 0.14, and at a triple point Xt — 3.15 where two ne- 
matic gels coexist with an isotropic gel. The coexistence of 
two nematic phases has been predicted in the literature in 
systems different from the one we have considered in this 
paper. In a solution of rods in the absence of an elastic 
network, it was predicted by Flory |3^ , |35| , and more re- 
cently it was predicted in Ref. pl[ in a system composed 
of a nematogenic network immersed in a nematogenic sol- 
vent. The coexistence of two nematic phases is a novel 
feature of our model, which is more subtle than in pre- 
vious studies of nematic gels since in our case only the 
network has nematogenic properties. The phase diagram 
of Fig. is shown in the coordinates osmotic pressure 77 
vs. x- As compared with Fig. ^|, the critical point is now 
on the other side of the line of coexistence of the isotropic 
and nematic phases. 



(a) 




X 

(b) 




X 

Fig. 7. Phase diagram for model 2 in the coordinates: (a) (j> vs. 
X and (b) 77 vs. \. The parameters are b/a = 100, f3 = 0.78, 
v c — 7 ■ 10~ 6 and <f>o = 3 • 10~ 3 . The collapse occurs at x ~ 4. 
The coexistence region of two nematic gels denoted N\ + N2 
terminates at a critical point located at x'c — 3.1 and 4>' c — 0.14. 
At Xt — 3.15 there is a triple point represented as a dotted 
line in (a) and by the point t in (b) where two nematic gels 
with different volume fractions and nematic order parameters 
coexist with an isotropic gel. 

4 Summary and Conclusion 

In this paper, we have studied the phase behavior and the 
collapse of a lyotropic nematic gel using Onsager theory to 
describe the isotropic-nematic transition. Upon decreas- 
ing the quality of the solvent sufficiently, we find that 
a nematic gel always undergoes a discontinuous volume 
change, even when the rods embedded in the gel are not 
physically linked to the network. The discontinuous vol- 
ume change is accompanied by an isotropic-nematic tran- 
sition, and for this reason it is of a different nature than 
the one discussed by Tanaka for isotropic gels. We have 
discussed the possible phase diagrams that these systems 
should exhibit; in particular, we have shown that three dis- 
tinct coexistence phases are possible, i.e. isotropic/isotropic 
(I/I), isotropic/nematic (I/N), and nematic /nematic (N/N). 
The I/N coexistence is always present, the I/I coexistence 
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and N/N coexistence may be present depending on the 
precise value of the parameters. 

Our model uses the Onsager free energy for rods in 
an isotropic solvent to describe rods in a polymer gel. We 
considered two models: Mmodel 1, in which rods are a part 
of the network and model 2, in which the rods are free to 
move within the free space of the gel but not permitted 
to escape from it. In model 2, we retained the full transla- 
tional entropy of the Onsager theory, whereas in model 1, 
we set it equal to zero since the rods formed a part of the 
polymer network. Thus, in both models, we overestimate 
the rotational entropy of the rods and we overestimate rod 
translational entropy in model 2. It would be interesting 
to consider models that provide more realistic treatments 
of these entropies. Our treatment also neglects completely 
the random torques and forces that the randomly cross- 
linked network exerts on the rods. These random fields are 
definitely present, and they may make it difficult to create 
an aligned lyotropic nematic gel in the laboratory. 

We hope, nevertheless, that this paper will motivate 
experimental work on gels with rods embedded in them. 
These experiments are important because of the poten- 
tial applications of these materials: these gels are in some 
sense a new type of sensors, which respond to a change 
in osmotic pressure by a change of volume. In particular, 
our work is relevant to some experiments which are now 
undertaken using crosslinked diblock copolymers as rigid 



rods embedded in a polymer network 37 
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